m tm i ^ s2 '° 


MAXIMUM POSSIBLE ERROR 
IN POSITION IN A LEAST 
SQUARE ORBIT 


GPO PRICE $ 
CFSTI PRICE(S) $ 


R. G 


BY 

. LANGEBARTEL 


$ N67-3Q145 

Z (ACCESSION NUMBER) 


(ACCESSION NUMBER) 

(thru) 

Jt/f' 

/ 

(PAGES) 

(CODE) 

_ TjM - 


(NASA CR OR TMX OR AD NUMBER) 

(category) 


Hard copy (HC) 
Microfiche (MF) 


ff 653 July B6 



2 *. && 


FEBRUARY 1967 


GODDARD SPACE FLIGHT CENTER 

, GREENBELT, MARYLAND 



MAXIMUM POSSIBLE ERROR IN POSITION 
IN A LEAST SQUARES ORBIT 

by 

R. G. Langebartel 


February, 1967 


Goddard Space Flight Center 
Greenbelt, Maryland 
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R. G. Langebartel 


SUMMARY 

A second order theory is developed for the problem of determining the 
maximum possible error in position if a Keplerian orbit is fitted by least squares 
to a set of observational data under the condition that the sum of the squares of 
the distances between the true and observed positions be held constant. A re- 
sult of the first order theory is that the maximum error at any one position in 
the orbit occurs when the observed positions coincide with the least squares 
computed positions. This doesn't remain true in the second order theory. 

Application of the second order theory is made in detail to the case of cir- 
cular orbits. 
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LIST OF SYMBOLS 

a semi-major axis of true orbit 

e eccentricity of true orbit 

P angle of perigee of true orbit 
z instant of perigee passage in true orbit 
a semi-major axis of least squares orbit 
e eccentricity of least squares orbit 
angle of perigee of least squares orbit 

£ instant of perigee passage in least squares orbit 
(r, f ) polar coordinates of position in true orbit 

(p, 9) polar coordinates of position in least squares orbit 

s mean anomaly for true orbit 
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<y mean anomaly for least squares orbit 
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/dj.b. see (1.4) 

W 2 see (2.1) .. 

Q 2 see (2.2) 

P see (3.28) 



MAXIMUM POSSIBLE ERROR IN POSITION 
IN A LEAST SQUARES ORBIT 

R. G. Langebartel 


INTRODUCTION 

A pertinent problem in orbit determination theory is the gauging of the effect 
of observational data errors. One approach, suggested by Dr. B. Kruger, is to 
investigate the nature of the distribution of errors that gives rise to the maxi- 
mum error in one predicted position. That is the problem considered here for 
the case of two-dimensional Keplerian motion where it is assumed the computed 
orbit is an ellipse fitted to the observational data by least squares. This least 
squares ellipse is, of course, a specialized curve fit in that it takes into account 
the special parameterization of the ellipse (with respect to time) peculiar to 
Newtonian two-body motion. As constraint in the maximum problem it is assumed 
that the sum of the squares of the distances between the observed and true posi- 
tions is a constant. No other assumption on the errors is made. 


THE LEAST SQUARES ORBIT 

The elements of the true elliptic orbit in the plane are the semi-major axis 
a , the eccentricity e, the angle of perigee p , and the instant of perigee passage 
z . The elements for the least squares orbit are indicated by the corresponding 
Greek letters: a. , e , n , and £ . Position in the true orbit is given by the polar 
coordinates r (distance) and f (angle), position in the least squares orbit by p 
and 9 . The direction from one focus as the origin from which angles are meas- 
ured is arbitrary but fixed throughout the discussion. The functional relations 

r r = r (a, e, p, z) fp = p (a, e , it, £) (1.1) 

\f = f (a,e,p,z) \cp = 9 (a, e ,77, £) 


1 


are those of the Newtonian theory: 


r 


< 




r = a (1 - e cos g) 
cos (f -p) = a r -1 (cos g - e) 
s = \/~iTa~ 3/2 (t - z) 
s = g - e sin g 


p = a (1 - e cos y) 
cos (cp - 7 r) - a p~ l (cos y - e). 


< 


(1.2) 

cr = \/Jl a~ 3/2 (t - £) 
cr = y - e sin y 


Thus, s and g are the mean and eccentric anomalies, respectively, for the true 
orbit, andcr andy are those for the least squares orbit. The gravitational para- 
meter p is the product of the gravitation constant by the sum of the masses of 
the two bodies. 


The least squares orbit is defined by requiring a, e, tt, and £ to be so 
chosen as to minimize the sum of the squares of the distances between the 
observed and computed positions. Suppose that positions have been observed 
at n + 1 different instants of time and let these observed positions be denoted 
by (r k ,f k ), k = 0, 1 n. Consequently, the function to be minimized is 


>2 = 


s ( ? k + Pi - 2 K Pv cos <^k - <Pk» ' 

k s0 


(1.3) 


Introduce the notation 




V. 


b 1 = a 

b 2 * e 
b 3 ^p 


b 4 = z 


r 


< 


/3 X = a 
^3 “ 71 

/3 4 = £ 


(1.4) 


The function P 2 is a function of the /3. through the variables P k and <p k . Hence 
the four equations of condition for the least squares orbit are 


1 

2 3 y 8 i 



(h-90] 


k 


f k p k sin (? k -9 k ) ^| = 0. (1.5) 


The solving of these equations determines /3 i in terms of the quantities r k and 
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2. THE MAXIMUM PROBLEM 


We wish to ascertain what distribution of errors in the observational data 
will bring about the greatest discrepancy between the true and least squares 
computed positions at the k th instant, t k , if we impose the constraint that the 
sum of the squares of the distances between the true and the observed positions, 
i.e. the sum of the squares of the errors, is held constant. That is to say, we 
are to maximize 


w2 3 r k +^k - 2r k^k cos ( f k -'Pk) ( 21 > 

subject to the side condition 

n 

Q 2 s YL { r k + ? k - 2 r k K cos ( f k - * k ) } = c2 • ( 2 - 2 > 

k=0 

The problem is to determine the values of r k and f k that make W 2 a maximum. 
The variables p k and 9 k in W 2 are to be regarded as functions of ? k and f k by 
virtue of (1.5) and (1.2). It should be noted that s, g, <r,y vary with t so we 
write 


r 


r k = a (1 - e cos g k ) 
cos (i k -p) = a r " 1 (cos g k - e) 

s k " ^a -372 (t k -z) 


s k = g k - e sin g k 


C 


< 


p k = a(l - e cos y k ) 
cos (tp k - 77 ) = ap -1 (cos "y k - e) 
cr k = /JIa ~ 3/2 (t k -.£) 


(2.3) 


cr k = y k - e sin y, 


k* 


The fundamental equations for the isoperimetric problem are 

dp. 3 m 

[^-r k cos( ?k -f k )] — 1 + r k/ 9 k sin (cp k - f k ) ~+ U. - r ra cos (f n - fj] =0 

3 r 3 r 


dp, d Op. 

[^k- r k c °s( 9 k - f k )] — +r p sift( 9 k -f k ) \ [r ra r in sin(f m _f m )] =0 

d f Bf 
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where k is the Lagrange multiplier. 


A more compact form for the equations involved results if we introduce 
the notation 


(K - r k cos (9 k - f k ) 

\\ = r k ^k sin (9 k - f k ) 

r k cos A - f k) 

\\ - r k ? k sin (f k - f k ) 

IV*-?k cos (9 k - K ) 

IA s f k ^ sin (9 k - K)- 

The necessary conditions for a maximum thus have the form 


(2.4) 


r 


3 9 k 

0, ^L + $ k __L +X0in = o > 

or dr 


m = 0,1, • • • , n. 


(2.5) 




3/° k ^9 k 

^-r + ^7r +k *« = 0 ’ 

Bf B f 


The dependence of p k and cp k on r m and f m has been specified by the least 
squares equations of condition (1.5). Use can be made of (1.5) if we write (2.5) 
in the form 


r 







= o 




V. 



^9 k ~ 


3/S. 

^ 9 m = 0 

B f 


(2.6) 


and then use (1.5) to eliminate 3/3./B r and d/3 /df . 

i m 'l' m 
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The equation (1.5) is 


V fn BPk * 39k 1 n 

Z_. h M + k 3/3. J" 


(2.7) 


k=0 


and upon differentiation with respect to r m and with respect to f m this becomes 

r 


V -1 °Pi 

i = l m 


j =1,2, 3, 4 


4 


i = l 


Wi 

* 

3f 


= K. 


jm 


where 


G ii" 


v - * r - 32 P\ 


d 2 f k 30, 3p, 34>, dcp k 


J J Cl K , *f\ K K'K K. K 

Lj Y* 3/3. 3/3. + k 3/3. 3/3. + 3£~ 3#" + 3£~ 3^ 


« 3p « 3<p 

H . 3 cos (f - <jp ) — + p sin (f - a>) ■■■■ ”- 

jm V m T m / ^ a 'm v m T m' 


(2.8) 


(2.9) 


3 /3 - 3 <jp 

K jm 5 - K Sln ( f m - 9*) + r m COS ( f m " 9 B ) ^g- 

L j j 

Upon introduction of the matrices 
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we have for (2.6) and (2.8) 

P k M + \ ® = 0 

B f 

{ ( 2 . 11 ) 

p M + K $ = o 
Bf 


gM = h 

3 r 


gM = k . 

Bf 


( 2 . 12 ) 


Evidently the necessary condition for a maximum now has the form 


fa G _1 H + K® = 0 
|P k G~ 1 K + X $ = 0 


(2.13) 


These equations are highly non-linear in the r m and f m and recourse must be 
had to some approximation procedure. Suppose that the functions involved are 
expandable in powers of a variable 8 : 
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P k =P< 1 > 8 +p(2) 82 + ... 

0 = 0 ( l > 8 + 0 ( 2 ) 8 2 + ... 
$ = $(1) S + $(2) e* + • • • 




G = G( & ) + G( 1 > 8 + G(2) ?2 + . . . 


H = H<°) + H( J ) 8 + H<2) 8 2 + • • • 


(2.14) 


K = K<°> + K( x ) 8 +K<2> 8 2 + ... 


^ \ = \«» + \( 1 ) g + \ca) 8 2 + .... 

Note from (2.4) and (2.10) that if the observed positions coincide with the true 
positions, P k ,0, and $ all vanish (which is, of course, in conflict with the 
constraint Q 2 = c 2 if c > 0) . If 8 is a parameter such that the observed con- 
figuration becomes the true configuration when 8 vanishes then the constant 
terms in the expansions of P k , 0 , and $ must all be zero as is indicated in 
(2.14) . In view of the parenthetical statement above we must also have c and 
8 vanishing simultaneously. Actually, we shall ultimately take 8 to be a mul- 
tiple of c but this specialization is not necessary now. It is only required at 
this point that the vanishing of 8 bring about the coincidence of the observed 
positions with the true positions. 

We need the expansion of the matrix G _1 . 


G" 1 =F<°) +F( X ) 8 +F<2) 8 2 + ... 

/. GG' 1 = I = (G<°> + G(2) 8 +G<2) 8 2 + ...) (F<°) +F< 1 ) 8 +F<2> S 2 + • • •). 
= G<°> F<°> + (G<°> F< J > + G< 1 ) F<°>) 8 + • • • . 


This shows that 


G~i = p(0> _ p(o> q( 1) p(0) g + 

-i 


F(°) = G(°) 


(2.15) 
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The substitution of (2.14) and (2.15) into (2.13) leads after some matrix algebra 
to the infinite system of equations 


p(i) f(°) h<°> + \<°) 0< l > = 0 
pa> p(°) k<°> + \(°> $o) = o 


(2.16) 


pa)F(°) - pa)F(°)G( 1 )F(°)H<°> + p( 2 )F< 0 >H(°> + \<°) 0 ( 2 > + \a)@a> = 0 


Fj»F(#)K(n - pa)FC0)G( 1 )F< 0) K<°> + P< 2 >F<°>K<°> + )^ 0 >$( 2 > + \(D$(D = q 


3. THE FIRST ORDER THEORY 

The particular character of the first order equations (2.16) rests on the 
manner in which the low order terms in the expansions of r , f , p , and 
are involved. Since S has been specialized to the extent that r k -*r k ,£ | k ~‘r k , 
f k f k , <p k ^k> ®-® £-»0 and k = 0, 1 , . . . , n, the expansions are written as 


; k = r k+ ?<‘>g t ;c*>e* t 


f„ = f k * f£ 2> e 2 + 


p k = r k+ p(‘>g^P>e 2 + --- 


<Pk = f k + <p£ 1> E + fk 2> ' 


Further, 


l£± (2) 02 

3 p 3b + L^J + L ^/5 J + 


1%. = H k 
3/3 3b 


M (1) g + M (2> 8 2 

L.3/3 J L^/SJ 
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where 




0^ = ©< 1 > e+0^2) 82 + ... 

% = «£*> e + <f^2) ga + .... 


(3.4) 


The first few terms in the expansions of the quantities appearing explicitly in 
(2.13) are easily obtainable and in particular we find 






iii 

Bb 


B f I B r J 
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Bb 


F<°) = 



r 2llll 

Bb Bb + * k Bb Bb 


-1 


H<°) = — 
Bb 





•Ha 

3b 4 / 
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K<°> = 




( 3 . 5 ) 


©(!) =X<1) - (?<»> , • • • , ?<D) 

$(D = x( j ) R S ( r 2 f(l) , ... , r 2 f(D) 


where 


^■('O H"' ••• ' r n H l) ) 


I T r 


R = 


‘n / 


There will also be needed the expansions of the least square elliptic orbit 
parameters. Upon writing 



/M 


/ b A 


r\ 

yS* 


b = 

• 

/30) ^ 

• 


W 1 

> 

\ b J 

> 

Wv 


we have 


/3 = b + /3( 1 ^ 8 + /3( 2 ) S 2 + • • • . 


( 3 . 6 ) 


( 3 . 7 ) 
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Now 

r 


fit 


(i) _ 


da 


dS 


J S = o 


yi ^a 

3/3. 

i*l 1 


3 a d^ 

d8 


K 1 ^ r U 

U1 1 


t 3 r^ 
3b 


(3.8) 


% 


(i) 


- M . ft ,12k d ^‘ 

L dE J8.o' a 3/5, 

1 si 


d8 


S*o 


\ 4 i 3 f, 

i si 1 


T 3 f. 


3b 


Consequently, there is the alternative expression for P£ 1 } : 


P(D =/3(D 


T / 3rl r2 

3b 3b + Fk 3b 3b 


(3.9) 


It is apparent that there is still needed the relation between ft 1 >and r ( 1 > , 
fW. This relation results from the lead off term in the expansion of (1.5), the 
equation of condition for the least squares orbit, in conjunction with (3.8). The 
coefficient of 8 in (1.5) when set to zero gives an equation which when rear- 
ranged has the form 


n 

L 

ksO 




3fT 




3b 


= £ r 


3 r l . 3 f J 

— — + rj* f^ 1 ) — 
3b k k 3b 


(3.10) 


And so the required relation between ft 1 > and r < 1 > , f ( 1 ) is afforded by 


ksO 


3 r. 3 r? 


3 f , 3 


+ r 






3b 3b k 3b 3b, 


< 311 > 


k.O 


The matrix R that appears in <I> (1) is an encumbrance that can be made to dis- 
appear from the equations (2.16) by the introduction of a new matrix J (0) in place 
of K (0) where 
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Note that K (0) = J (0) Rwith R non- singular, so that the equations (2.16) have 
the form 


p(D F (o> H <o) +x( o) X (i) = 0 


p(l> F (0) j(0) + X (0) X (i) = 0< 


It may be noted that B r k / B b and r k B f k /Bb are the k + 1 st columns of the 
matrices H (0) and J (0) , respectively, whereupon there is suggested the notation 



And it is convenient to combine these into one matrix: 



(3.14) 
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The matrix of these matrices likewise occurs and so we define 


_BX _ 

Bb 


Then 


G(°) 


/ Br ° 


B r i 


/ ^ b l 

• 

BW 

r 

1 Bb x 

1 3r o 

• 

* f 0 


r Bf i 

Qj 1 

O - 

1 

r ° ^b 4 

3b 4 

Bb 4 

n 




= £ 

(H^°> HC°) T 

+ J< 

°> J^ 0)T ) = 

TIo 




■up 

ii 

si* 

II 

BX 

Bb 

BX T 

Bb ’ 

ii 

'Co 

/— s 

I-* 

H 

BX k B Xl 

Bb Bb ' 




Br 


Bf 


Bb. n Bb. 


B r 


Bf 


• • * 


Bb. n Bb. 


/3(D t = (xjO H^ 0)T + XC 1 ) J(°> T ) F<°> 


(3.15) 


(3.16) 


(3.17) 

(3.18) 


Thus the observed position coordinates as represented by the matrices and 
X^ 1 ) must satisfy 

3 X 3x t 

CXd) h(°) T + X'b t( 0) t \ F (0) — Ji — i F (0) H (0) + \(0) x(D = 0 
^ i " j > Bb Bb 1 


< 


(x<‘> h<°> t tX^) J(°) T )F(° 


3 X 3 yT 

) * i F (0) j(0) + \(0) X(1) - 0 


(3.19) 


Bb Bb 


This is an eigenvalue problem for the eigenvectors X£ x) , XS, 1 ) and the 
eigenvalue A. (0) . Some information about the eigenvectors is readily obtained by 
postmultiplying the first equation by H (0)T , the second by J (0)T , and subtracting 
since these operations, if f 0 , result in 
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(3.20) 


X(D h<°> T (H<°> H( 0 ) T ) _1 - J<°) T (J<°> J (0)T ) _1 = 0. 

This cannot be used by itself to solve for X® in terms of X^ 1 ) or vice versa if 
n > 3 because H (0) andj<°) are non- square and thus have no inverses. 

The two matrix equations (3.19) represent 2n + 2 equations in the unknowns 
r (! ) , f O ) and so considerable complexity in the spectrum could be expected. 
Fortunately, this is not the case as only two of the eigenvalues can be different 
from zero. This circumstance arises essentially from the fact that BX k /Bb is 
a non-square matrix. Let X s (X^ 1 ) ,X^ 1) ), i.e. a row matrix with 2n +2 elements, 
and let L be the 2n + 2 x 4 matrix 


Then the matrix A = LF (0) BX^Bb is a 2n +2 x 2 matrix and the eigenvalue 
equation becomes simply 


xaa t + x<°> X = 0. 


The two columns of A may be regarded as the components of two vectors p 
and q in a vector space of 2n + 2 dimensions. The vectors orthogonal to both 
P and q fill out a 2n -dimensional vector space. Let X 0 be any non-zero 
vector in this 2n -dimensional space. Then X Q AA T = 0A T = 0 and so every 
vector in this 2n -dimensional subspace is annihilated by AA T and is therefore 
an eigenvector of this matrix with eigenvalue = 0 . Since the 2n -dimensional 
space is spanned by 2n independent vectors the multiplicity of the zero eigen- 
value is 2n . In the unlikely event that p and q are not independent but are non- 
zero the vectors orthogonal to bothp and q will fill out a space of 2n + l dimen- 
sions and the multiplicity of the zero eigenvalue in this case will be 2n + l . 

Thus, in general, there are only two non-zero eigenvalues. The search for 
these two eigenvalues is considerably simplified by reducing the eigenvalue 
equation in X ^ , X® to one in/l (1) . This is accomplished by postmultiplying the 
first equation in (3.19) by H (0)T , the ’’ wfV ' T 
(3.18) and (3.16). This eigenvalue 

/BX, BXJ 

\ Bb Bb 


second oy j , ana tnen adding ana using 
equation for /3 (1) , 


:o> ) 


+ \(°> G<°>) ySd> - 0, 


(3.21) 
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can be analyzed in the manner above to obtain the result that in this fourth 
order system two of the eigenvalues are, in general, zero. The zero eigenvalues 
correspond to the eigenvectors that are orthogonal to botn columns of 
3X k /3b,i.e., 3 r£ /3b = r fe 3 fj /Bb /3 ( - 1 '>= 0. But this implies by (3.8) 

that = r k Tk 15 = 0 which in turn implies, as we shall see later, that at least 
to the first and second order W 2 = 0 , indicating, a .minimum jrather than a maximum. - 

The eigenvalues are, of course, the roots of the determinantal equation 


US. . x,o) m MI 

Bb Bb Bb Bb 


= 0 . 


(3.22) 


The eigenvalues are all real since 3X k /3b 3x£ /Bb and BX/Bb BxV^bare 
both symmetric with real elements [ref. 1, p. 306 ] . 

We now show that the original eigenvectors, X^ , X^ , have a very simple re- 
lationship to tne p< l \ cp^ . If (3.20) is used in (3.18) there is obtained 


/3(1) T = [xjn H<°) T + X< x > H<°> T (H<°) H<°> T ) -1 (J< 0) J (0)T )] f<°) 
= X^ 1 ) H^°^ T (H^°) H^°^ T ) _1 [H<°> H^°^ T + J<°) J C 0 ] p(°) 

= x^ 1 ) h(°) t (h^°) h ^ 0)T ) _1 ^ 0) ^ (0) 

= x [ 1 > h(°) t (h<°> h<°> t ) -1 . 


Similarly, 


Hence, 


/3(D T =X(D J(°> T (j<°> J(°> T 


"x(i) H(°> = /3(») H<°> H<°> , 

X(i) j(°) T = p(i) r J(0) J(0)T 


(3.24) 


(3.25) 


We cannot conclude from this alone that X^ x) = /^ X) H (0) and X^ 1:> = /3^ 1) J (0) 
since H (0) and J (0) , being non-square for n > 3 , have no inverses. However, 
these formulae are indeed true for eigenvectors corresponding to non-zero 
eigenvalues. We can conclude from (3.25) that X^ x) is H (0) except possibly 

for an additive zero divisor of H(°> T . That is, we can write 
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X i (i) = /? (i) H <®> f Y l , 


x(D =/ e(D J(0) + y 2> 


where 


Y x H° = 0, 


Y 2 J° = 0 . 


Substitute these into the eigenvalue equations (3.19) to give 


r 




d) T ( F ( o ) x(°) lW 0) + X<°> Y, = 0, 

\9b 3b J 1 


< 


/3 ( 


i) T l^L ^f<°) + \<°) 1^ J(°) + \(°) 
\ Bb 2h J 


y 2 = o. 


But if X( J ), X^ 1 > are eigenvectors then by (3.21) 

^ (1) T p(0) x(°> i) =0. 

P \ Bb 3b ) 

Therefore, if 0, we must have Y t = Y 2 = 0 and, consequently. 


X(D = /?(i) T H(°>, 
X(i> =/3 (i) T J(0). 


(3.26) 


This joins up with (3.8) to give the result that if A.(°>^0, 



(3.27) 


that is to say, for a maximum to occur the observed positions must coincide with 
the least squares computed positions and therefore must fall along a Keplerian 
ellipse (though not, in general, the true orbit). This property, to be sure, is a 
consequence of the first order theory. As we shall discover later, it does not 
carry over to the second order theory. 
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Nothing has been said so far about the normalization of the eigenvectors. The 
size of the eigenvectors is regulated by the side condition Q 2 = c 2 . Before we 
can determine its effect we must decide on the variation of c with S . We choose 


8 = 


a /n + 1 

Since the lead off term in the expansion of Q 2 is 


_ - - -< 3 , 28 > 


£ ( f k‘ >! + r ? ? £ 1>2 ) 82 ■ 

k=0 ' 


the eigenvectors must satisfy 

X( J ) Xp) T + X( 1 )X 2 (D T = a 2 (n + 1). (3.29) 

The vector ft 1 Hs governed accordingly. Working with (3.26) we find 

X< 1 > X^^ + XO) XO) T = /3(i) T H<°) [ft^ H(°)) T +/ 5( 1 > T J(°) [ft 1 ^ j(°)j T 

- ft l ) T H<°) H(°) T ft 1 ) + ft 1 ^ J(°> j(°) T ft 1 ') 

= /3C) T G<°> ft l \ 

Consequently, /3Wis normalized by 


fit? fit) = *2 (n + 1) . (3>30) 


The maximum squared distance between the k + 1 st true point and the k + 1 st 
least squares computed point is 

W " * r k + Pi ~ 2r k P k cos -V 

(/?<d T H<»)) T + (/3<‘) T J<°>)(/3“> T j{ 0 >Vle j 
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= /?(i) T (h(°> H(°) T + J«» J(°) T ) /3(D 8 2 + • • 


t bx u 3X. 

= /3(D — i ! ySd) S 2 + . . . . 


Bb Bb 


(3.31) 


This can be expressed in terms of the eigenvalue \<°>by seeing from (3.21) that 


T BX U 'dXj T 

/ 3 O) — Jl — L £(1) _ _ \( 0 ) ^(i) G ( 0 ) ySfO 

Bb Bb 


and this in turn is -\<°> a 2 (n + 1 ) by (3.30). Hence, an alternative form for 


W 2 is 


W 2 = - \<°> c 2 + . . . 


(3.32) 


4. THE SECOND ORDER THEORY 

Pushing the analysis to the second order consists essentially in treating the 
set of equations (2.17). The following matrices will be needed in the work: 


X< 2 > - 


• , ?< 2 >) 

n 


X<2) = /- r f(2) 

2 ' r 0 0 ’ 


, r n 


y(i) = 


/ 




f(i> 

0 


f(l) 

1 0 


f(Dy 

n 


r (i) 

0 


r(D = 


f(i) 

1 


\ 

0 

';<■>/ 


(4.1) 
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Then 


r§(2> = x(2) + lx(D yd) 

[$( 2 )r- 1 = X< 2 > +X[V v< x > 


(4.2) 


'-"■[nr 

■•”*■( [3] 


(1) R + 2 RW] r 

ob j 


G (1) = H d) H (0 ) t + jd) J(0) T + || [|| 

J p(2) = Ad) + P p H(°) T + r k 9 p J(O) 1 


H Id) 1 


(4.3) 


A £ 1} - 


V. 


pivl 0 ’ 1 1 m 

LsjJ +r 2 ^L 5 jJ 


a 3 r J a II 

* * n\ _v ia k 


Bb 


+ r k 


Also, from (3.17) and (3.21) we get 

P (i) _ _ x(0) ^(D 1 g(°). 

The fundamental second order equations (2.17) thus have the form 


(4.4) 


- A<°> /3d) T Hd) + A<°> / 3d) T G< J > F<°> H<°> + Ad> F<°> H<°> + (H< k °> T p< k 2 > + J< 0)T r k <?< 2 )) F<°> H<°> 
+ \(°) ^X<2) + Ix<D yd) j + Ad) ySd) T H(°> = 0 

- A<°> /3d) T j(D + \(0) /3(D T Gd) F<°> J (0) + Ad> F (0) J(°> + (H<°) T pp +]p T r k cp< 2 >) F<°) J<°) 

_+ \(°) (X(2) + xd) Vd) ) + Ad) /3 (D t jW = 0 (4.5) 


These equations involve two sets of second order unknowns, p ^ 2 ^, and X[ 2 \ 
X< 2) , besides the eigenvalue perturbation Ad), in order to obtain the necessary 
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further relations among these we note first that 


! 


■s[ 


d2 ^k " 

dS 2 - 


8 = o 


1 

2 


d y» d P k d A 
dS l_j 3/T dT 

i = 1 1 


£*o 


3r, 1 v^n ir—\ B 2 r, 

-L ^ 2>+ l 

i=l 1 i=l j=l 1 1 


with a similar formula for <p£ 2) , and that 


3 A _9r 
3/3. 3b 


7*L£k*>*-~ 

1 W i 1 


so that 



3 2 r v 

— yS? 1 ) 

3b. 3b. J 



The connection between /o£ 2) , 9^ 2 ) and the orbital elements is thereupon seen to 
be, in matrix form, 


Pi 2) =/ 3 ( 2 )T h£°> 


V. 


r k^ 2) 


= /S (2> T J(0 } 



(4.6) 


To complete the connection of pC 2 ),^ 2 ) with x£ 2) , X^ 2) we need now the relation 
between /3< 2 ) andX^ 2) , X^ 2) . This relation resides in the least squares equation 
of condition (1.5) which, when expanded to the second order, gives 


/3< 2 ) T G<°> = X< 2 > HC°> T + X< 2 > J(°) T + SI, 



BSld) 3 XT 

3/S J 3b 


(4.7) 
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the counterpart to (3.18). We note further from (4.6) that 


T T T 3X. BX? 

P ( 2 )H(°) + r k ^)j(0) =/3 (2) ^ ^-a k 


(4.8) 


n^-I/^) 7 




w 


3b 


This enables us to remove p®\ ? k 2) from (4.5) in favor of /3< 2 >, thereby giving 


C t 3X. BX? 

bw — i — i. fw h<°) + a<°) x$ 2 > 

^ 3b 3b 1 


= A<°> /3d) T h<» - A<°) /3d) 7 Gd) F<°) H<°) - Ad) F< 0) H<°) +G k F<°> H<°) 

_ a<°) lx ?) v d> - Ad) /3d) 7 h<°) * n , 

2 2 1 

(4.9) 

T BX k BX* 

ySC 2) i — ii F<°) J<°) + \<0> x< 2 > 

Bb Bb 2 

= \(°) /3d) T jd> - \(°) /3d) T Gd) F (0) J<°> - A< J ) F<°> J<°> + fl k F<°> J<°> 

- a<°) x w v - A( J ) /3d) 7 j<°) ^ n 2 


The unknowns Xd), X< 2 ), can be eliminated from these equations to leave an equa- 
tion in /3d) alone by postmultiplying these two equations by H ( °) T and J (0)T , re- 


spectively, adding, and using (4.7). The resulting equation is 


/3(2)‘ 


,3b 3b 


+ \<°)g<°)) =n 1 H(°) 7 +n 2 j(°) T + A(°)n, (4 .io) 


the counterpart to (3.21) . Since Ad) has already been chosen so as to make the 
determinant of the coefficient of /3d) T vanish, the number Ad) , which is present 
in n x and II 2 , must be chosen so as to make the rank of the augmented 
matrix equal to that of the coefficient matrix. After this has been done, (4.10) 
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can be solved for /3 (2) . If, on the other hand, /3< 2) is removed from (4.9) by the 
use of (4.7) there are displayed equations from which x[ 2) , X^ 2) can be obtained. 

A simple relation exists between X| 2) and X< 2 ), a relation derivable by 
multiplying the first equation of (4.9) by H(°) T (H<°) H(°) t ) 1 , the second by 
j(°) T (j(°) J (°) T ) -1 , and subtracting. It is 

x< 2 > h(°) t (h<°> n^ r y' - x(2> j(°) T (j<°> jo) 1 )- 1 

= yS(l) T _Ij(°)vW)H<°> T (HWhW 7 ) -1 - (jd) -H(°> VW) JW T (J(°) j(°) T )~l] . 

(4.11) 


This is evidently the counterpart of (3.20). 

The coefficient of £ 3 in the expansion representing the constraining equation 
(2.2) is 


2 


r. ? ^ 2) + r k vp +7 

kzO 



so that we must have 


X^X^' + X^X^‘ +Ix< 1 > (X< 1 )V< 1 )) T = 0. (4. 12 ) 


which is for the second order theory what (3.29) is for the first. The cor- 
responding equation for /3d) results from the application of (3.26) and (4.7) to 
(4.12). We find that 


/3< 2 > 


T ML MI /?( i) + i bw t ML 

3b 3b H 2 3b 


' W /3d) + 1 /3d> T JW Vd> HW T /3d) = 0 , ( 4.i3) 


the counterpart to (3.30). 

A similar computation leads to the expression for W 2 , 
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w 2 = ( p <!) 2 + r 2 cpCl> 2 ) S 2 + (2^> 4 2 > + 2 r 2 ^ <p< 2 > + r k cp(^) £3 + 

b/ kd t ^ !S^) 8 . 

^ 3b 3b 


+ 12 


T BX b BXj T BX. T BE*. 1 (1)T 

V ^[w\ * 


i) 


+ /3(’> t j<°> f< x >H( 0 > pmy 8 3 + .... 


(4.14) 


Note that the coefficient of £ 3 (as well as that of £ 2 ) vanishes if /3 (1 ) T BX k /Bb = 0 
which occurs for \( 0 ) =0. Thus the vanishing of the eigenvalue causes the vanish- 
ing of W 2 through terms in £ 3 at least. 


5. KEPLERIAN FORMULAE 

The work of the last two sections is independent of the type of orbit involved. 
All that is required is that the orbit depend in a well-behaved fashion on four 
parameters (the extension to any number of parameters is immediate). In the 
case of two-dimensional Newtonian two-body motion the parameter dependence 
is given by (1.2). From these the fundamental quantities appearing in the for- 
mulae of the preceding sections are easily computed. The column matrices 
h£ 0) and j(°> are 


A" 


1 r k "| a r k _1 e s k sin e k \ 


H<°) = 


- a cos (f k - p) 
0 


- /JI a 1/2 r " 1 e sin g k J 


(5.1) 
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- a r 


/ttp 


j^°> = 


[a + r k (1 - e 2 ) -1 ] sin (f - p) 


- yjl a 1,/2 r~ 1 / 1 - 


The matrices H (0) and J<°> are 4 x n + 1 and have the elements of H£°> 
J(°), respectively, as their k + 1 columns. 

Let 0) represent the element in the i th row and j th column of the 
matrix G(°). Then 


n 

= H [ a ~ 2 r k - 3 e s k sin gk + J a r k X (2 - a_1 r k) s k]> 

k = 0 


G iv = £ [ - 3a2 ^ 1 Sk sin gk " Fk cos (fk ' p)3 ’ 

k = 0 


<*s> 


= [-§./rn* s k ], 

k = 0 


n r 
* I 


^ |-//I a _1/2 e sin g k + a 1/2 r" 1 (2 - a -1 r fc ) s fe j 

k=0 


G<°> 
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n 

= [a 2 + a 2 {2ar“ x + (1 - e 2 ) -1 } sin 2 g fe ] , 

k=0 


_ L [ar k + r 2 (1 - e 2 ) -1 ] sin (f k - p). 
k=0 


(5.2) 


and 


4x4 
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n 



k=0 


n 



k = 0 


L [-✓£ a 1 / 2 /FTP] , 

k=0 


G^ = Ua' 1 r^ 1 ( 2a ~ r k^ • 

k=0 


( 5 . 3 ) 


This determines all the elements of G<°) since G (0 > is symmetric. 

If we suppose that the observations are made at equal intervals in time then 
the summations with respect to k that are present in the formulae for G< 0 >can 
be replaced by summations of a different sort. This is accomplished by 
employing the Fourier series expansions in s for the various orbital quanti- 
ties, since then the original summation index k will occur only on the variable 
s and will appear only as a simple factor. The sum over k is then accom- 
plished through the use of the following series which represent the complex 
geometric series and its first two derivatives: 


r 


k=0 


- n + 1 


< 


k*0 


(n + 1) 


L*- l v 


(n + 1) (2n + 1) 


V. k=0 
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1 1 

/ cos pk 6 = esc — p 6 sin n + ^ pd cos — p 8 

Z__I 2 2 2 


n 

y k cos pk<9 = n esc - p d sin ^ + * p0 cos - p(9 

/ j 2 2 2 


- 1 esc 2 — p<9 [(n + 1) - n cos pd - cos np0] 

4 2 


n 

^ k 2 cos pk0 - n 2 esc i p0 sin ^ ^ p/9 cos p0 


- — esc 2 — p(9 [(n + 1) - n cos pd - cos np/9l 

4 2 


-esc 3 — pd n sin i pd • (1 + cos npd) - sin np d cos ~p/9 
4 2 L 2 2 J 


n 

y ' sin pk0 = esc ^ pd sin — ■ * * pd sin p0 


n 

/ k sin pk# = n esc i- pd sin — I pd sin — pd 

Z__i 2 2 2 


- 4 esc 2 i p# [n sin - sin np<9] 

4 2 


n 

y k 2 sin pk# = n 2 esc i pd sin ^-t^.pd sin ^-pd 

L i 2 2 2 


- jCSC 2 -jp# t n s ^ n P^ “ S2Ln n P^ 


i esc 3 i p0 j^n sin i pd sinnp6> + cosi pd -(cos npd - 1)J. 
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Let ^ be one of the quantities connected with the orbit. Its Fourier expansion 
in terms of the mean anomaly s is of the form 


*k= A + 



(a p cos ps k + b p sin ps k ) 


(5.5) 


where s k is the value of s at the k +lst instant: s k = ^ a 372 (t k - z ). The con- 
stants A, a p , b are all independent of k. We are assuming the observations 
are made at equal intervals in time and so we write 

s k = k0 - 0 (5.6) 

where 6 = i/JlZ a -3/2 r, 0 = /^a _3/2 z , and r is the interval (in time) between 
successive observations. Then the Fourier expansion of becomes 


X k =A + 



|[a p cos p© - b p sin p0] cos pk 6 


+ [a p sin p0 + b p cos p0] 


sin 



(5.7) 


For the elements of G ( °) we require the sum of ^ over k together with its 
first two moments. In forming this finite sum we add the Fourier series term 
by term, a process justified bv the fact that the Fourier series converge for all 
values of s. [ ref. 2, p. 210 ] . The result is, upon using (5.4), 


u = n 


= (n + 1) 


A + 


csclptf sin (a p cos ps n/2 +b p sin ps n/2 ) 

p=l 


where s n/2 is the value of the mean anomaly at the instant t=(n/2)r. In a 
similar fashion we can compute the first two moments. The formulae are 
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11 

( a p cos P S n /2 + S ^ n P S n/2^* 


k=0 


P *1 


n n ” 

r. "A = S n /2 22 \ + ^ £ Vp' (a p Sin PS n/2 - b P COS PS „/ 2 >> 

k=0 k=0 p=l 




(5.8) 


2>- 

k=0 


(s 2 - 2s 2 ) 


n/2 


n n 

*k + 2S n/2 ^ 

k=0 k=0 


V. 

where 


+ 02 ' 22 D P 2> ’ (3p C ° S PS " /2 + b p sin ps n/2^ 

p=l 


D^ 0 ) = esc — p 6 sin — + - pd , 


D^ 1 ) = — csc 2 ip0 - (n sin 1 p0 cos — + - p 6 - sin — pd Y 
p 9 2 V 2 2 2 / 

V ' (5.9) 

=icsc 3 ip6 1, n 2 sin 2 — pg sin^tilpff + n sin — pO • (cos — pd cos - + 

p 6 2 L 2 2 2 V 2 2 k 

+ 2cos^-p#j -3 cos ip0 sin^-p^J, 


r 


n/2 





s 2 = _ 'ifo, tH e 2 - ne © + 0 2 . 
6 


(5.10) 


In order to use these formulae to obtain representations for the elements 
of G (0) we need to know A, a , b , the coefficients in the Fourier expansions 

P P 

for the various quantities connected with the true orbit. These expansions are 
known [ref. 2, p. 205] and the particular ones we require are given in the 
following table: 
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sin g 

0 

0 

— J p (Pe) 
ep p 

cos g 

1 

e 

2 


0 

r 

a H e2 ) 

-2^j; (pe) 

0 

r" 1 

a~ l 

2a’ 1 J p (pe) 

0 

r 2 

a2 ( 1+ l e2 ) 

4a 2 T 

- — J P (pe) 

p 2 

0 

cos (f - p) 

- e 

2(1 - e 2 ) , N 

e J P (P ° 

0 

sin (f - p) 

0 

0 

2 /I - e 2 j; (pe) 

r cos (f - p) 

3 

' 2 36 

T r > (pe) 

0 

r sin (f - p) 

0 

0 

2a /l - e 2 . . 

J B (pe) 

ep p 

sin 2 g 

1 

2 

-4P - r> J p (pe)l 

d(pe) [pe p J 

0 

r“ l sin g 

0 

0 

2 a -1 j; (pe) 

r 2 sin (f - p) 

0 

0 

a 2 /l - e 2 ^ 2 e + p 2e J p (pe) + J_ J' p (pe)] 

r" 1 sin 2 g 

la” 1 

2 

2a ‘ J » (pe)+ oe (pe) ] 

0 


In these formulae J (z) is the Bessel function of the first kind and pth order, 
and J' (z) is its derivative with respect to the argument, z. Then, as alter- 
natives to (5.3) we have 

00 

G<°> = (n + 1) ^1 +| e 2 +-|s 2 ^ - 4 ^ D £ 0) P’ 2 Jp (P e > cos P s n /2 


00 00 

“ 6S n/2 L D<°>p _1 J p (pe) sinps n/2 + 9 s 2 ^ D £ 0) J p Cpe) cos p s n / 2 

p“l p = » 
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00 <P 

+ 66 L D(») p" 1 J p (pe)cosps n/2 + 18s n/2 0 D p‘ } J p (pe). sin P s n/2 

P=1 P=1 


- + 9 ^ D< 2 > J p (pe) cos ps n/2 , 

P=1 


00 00 

Gj 2 ^ = | (n + 1) ae - ^ D p 0) 2a P _1 Jp (P e ^> cos P s n /2 “ 6 as n/2 22 D p 0> Jp^ pe) sinps n/2 

P=1 P = 1 


+ 6a6 D^ 1 > J p (pe) cos ps n/2 , 

p=i 


G [ ® ) = - — a A - e 2 

13 2 


(n + 1) s n/2 , 


G<°> =|^ a" 1 / 2 (n + 1) s„ /2 - 2 * L D^°) p " 1 J p (pe) sin ps n/2 

p-i 


00 UU 

6 //Li a" 1/2 s n/2 22 D< 0) J p (pe) cos ps n/2 + 6 //x a _1/2 6 D^ 1) J p (pe) sin PS n/2 , 

P =i P =i 


=(n+l)a 2 


5 - 4e 2 
2(1 -e 2 ) 


+ 4a 2 Y' D<°> 


P“ 1 


-P(l-e 2) lli J (pe) + 1-62 V (pe) 


e 2 (1 - e 2 ) p 


e(l - e 2 ) p 


cos ps 


n/ 2 * 


G(o ) 


4a 2 

A - e 2 


22 D(°> p' 2 (pe) sin ps n/2 , 

p=l 


G 24'* = ~ 4 //X a 1/2 22 D p 0> Jp (P e ) sin P S n / 2 ’ 

P=1 


(5.11) 
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= (n + 1) a 2 ^1 +-| e2 j -4a 2 D<°> p ' 2 J p (pe) cos PS n/2 , 

p=l 


634 ^ = - v/Ju a 1/2 /l - e 2 (n + 1 ) , 

00 

G^ = (n + 1) fj. a -1 + 4fia~ 1 D<°> J p (pe) cos ps n/2 - 

p=i 


The factors J p (pe) and J' (pe) in the summands in these series bring 
about quite rapid convergence. In fact, these factors are both positive decreas- 
ing functions of p [ ref. 3, p. 254] and, moreover, 



(pe) ^ (277p)“ 1/2 (1 - e 2 ) -1/4 exp 
(pe) - ( 277 pe 2 ) _1/2 (1 + e 2 ) 1/4 exp 





(5.12) 


which shows that the convergence is exponential inasmuch as the coefficient of 
p in the argument of the exponential functions in (5.12) is negative for 0 < e < 1 
If e is very small the first term in the series for G( 0) serves as a good 
approximation for the series. Besides this only a small extra error will be 
introduced if Jj(e) and J' (e) are replaced by (l/2)e and 1/2, respectively, 
since the error terms are 0 (e 3 ) and 0 (e 2 ) , respectively. 

The quantities occurring in the second order theory for the maximum 
problem depend to a great extent on the second derivatives of r and f with 
respect to the orbital parameters. In a straightforward computation these are 
found to be 


— r = — -- = — a -1 e(l - e 2 )“ 1/2 s sin(f -p) + — ar ’ 2 e s 2 cos (f - p), 
Bb 2 da 2 4 4 


= - cos (f - p) - £ a 2 r -2 s /l - e 2 sin (f - p), 

ob. db, de da 2 
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B 2 r 

3b 3 % t 


B 2 r 
Bp Ba 


= 0, 


B 2 r 

Bb Bb, 

4 1 


4 ^-.5 ' - = i Sjl a * 3/2 e(l - e 2 )" 1/2 «in (f - P) + — V 
Bz oa 2 2 


^-1 = -L- !- = [a 2 r" 1 + a(l - e 2 ) -1 ] sin 2 (f - p), 
Bb 2 Be 2 


B 2 r 

Bb 3 Bb 2 


B 2 r 
Bp Be 


= 0, 


B 2 r 

Bb 4 Bb 2 


B 2 r 
Bz Be 


= - yji a 3/2 r“ 2 /l - e 2 


sin (f - p), 



B 2 r 
Bp 2 


= 0, 


B 2 r _ B 2 r 

~Bb 4 Bb 3 Bz Bp 


= 0, 



e cos (f - p), 


B 2 f 15 _ 2 n 2 9 2 . /f 

= — r 2 s v 1 - e 2 - — ar 3 es 2 sin (f - p) , 

Bb 2 4 2 


M a 1/2 r -2 es cos (f -p), 
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a2f 3 — a r~ 3 

3b 2 3b, 2 


s (1 - e 2 ) _1/2 


[re - 2a (1 - e 2 ) cos (f - p) ] , 


B 2 f 

Bb 3 Bbj 


= 0 , 


3 2 f 

3b Bb, 

4 1 


- — i/a a 

2 


r 


-3 


[ r /]T- e 2 - 2aes sin (f - p)] , 


^ - {a 2 r -2 cos (f -p) + [a r _1 + (1 - e 2 ) -1 ] 2 cos(f-p)+2e(l- e 2 ) -2 } sin(f - p), 
Bb 2 C J 


B 2 f 

Bb 3 Bb 2 


= o, 


B 2 f 

Bb, Bb 2 


= yji a 1/2 


r" 3 (1 _ e 2 )" 1/2 


[re - 2a(l - e 2 ) cos (f - p) ] , 


B 2 f 


= 0, 


B 2 f 
Bb Bb 

4 3 


= o, 


- - 2 fj. r -3 e s in ( f - p) . 

Bb 2 
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In the work leading to the expressions (4.6) for p (2) andr<jp( 2 >there are the 
relations 


(5.13) 


showing the connection of the second derivatives of r and f with the matrices 
occurring in section 4. Indeed, with /3 (1 determined in a particular case from 
the first order theory and with [3p/3/3.] (Oand [3<p/3/3.] (1 ) computed with 
the aid of (5.13) and the above second order derivative formulae all matrices 
are computable that are needed in (4.10) so that this equation may be solved 
for ft - 2 ) and >S 1 > . 


< 


V 

B / 3 

L 'u 


3cp 


-l<»> 


( 1 ) 


= E 


1 = 1 


J 1 


* 

■E 

j=i 


3 2 f 


3b. 3b. J 

j i 




6. CIRCULAR ORBIT 

In the case of a circular orbit the eccentricity is zero and the angle of perigee 
and instant of perigee passage become indeterminate. Thus, in dealing with cir- 
cular orbits we have only to consider the radius, a . Consequently, such matrices 
as °) and G^ 0 ) which are4x (n + l)and4x 4 in the two-dimensional elliptic orbit 
case become lx (n + 1) and 1 x 1 in the two-dimensional circular orbit case. 

The true circular orbit is defined by 


r 

f 


= a 

= //x a -3/2 t 


(6.1) 


and it is on these formulae that the matrices occurring in the maximum problem 
are based. In particular, 
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^ H(°) = (1, 1, 1, • • • , 1) 

J(°) = -2/^ a _3/2 r (0, 1, 2 n) 




sh (*• -f‘ /?a ' 3/2kT ) 

|2= U, O, 1, - 2 /jH a~ 3/2 r, • • • , 1, - 2 v/^6 a" 3/2 nrj 


G<°) = (n + 1) 




1 + 2 jx a -3 n(2n + l - ) r 2 J 


( 6 . 2 ) 


The position at the k + 1 st instant in the least squares circular orbit is 
determined by 


(1) = K(°> ad) = ad) 


^d) =H ( 0 )a 




( 1 ) = a/ 1 ) a ® ^ ki 


(6.3) 


which were computed from (3.8) and (6.2). The first order Langrange multiplier, 
\<°>, is found from (3.22) and (6.2) to be 


• (6.4) 

but the one-dimensional eigenvector /3d) = a < l )is not determined by the equation 
(3.21) precisely because the matrices involved are lxl. That is to say, the 
only freedom for the direction of / S< 1 ) is that afforded by the sign of /3 (1 > and this 
is left unspecified by (3.21) . The magnitude of fiV) is, however, another matter 
and it is found from (3.30): 


A/°) = - (l + 2 fj. a -3 k 2 t 2 \ (n + l)” 1 


1 +^- / aa -3 n(2n + 1) r 2 
8 


fl 1 ) - a( 1 ^ = ± a 


[-1 




~ 3 n(2n + 1) t 2 


• 1/2 


Probably the most satisfactory way of representing the maximum value of 
W 2 is in terms of the true anomaly, f . By (3.32), (6.4), and (6.1) we obtain for 
the maximum value of the squared distance between the k + 1st true point and the 
k + 1st least squares computed point (to the first order) 
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w 2 = 



3 2n + 1 
+ 8 n 



c 2 /(n + 1) . 


(6.5) 


To get an idea of the size of this we note that if k = n and f n = 1 radian with n 
regarded as large W 2 is approximately 2e 2 /(nj. 1). . - - 

In extending the work to include the second order terms we need the follow- 
ing basic results: 


r . -- -i. 




KD 


B/SJ 


= 0 




LB/3J 


(i) 


- //x a” 7/2 a/ 1 ) kr 


Hd) = (0, 0, 0, . ■ . , 0) = 0-H<°> 


j(i) r: /J1 a^ 1 ) a” 5/2 r (0, 1, 2, • . . , n) = a’ 1 


(6.6) 


= _^.n(n + 1) (2n + 1) /xa~ 4 r 2 
8 

A<»> = - — fj.a ( 1 ^ a -4 n 2 r 2 

n 4 * 


The second order terms in the least squares computed position at the k + 1 st 
instant are then available from (4.6) as expressions in terms of a/ 2 > and these 
are 


2 > = a< 2 >, 

<jp£ 2 ^ = VJl a“ 7/2 kr ) 



(6.7) 


38 


_ The remaining second order quantities present no problems and we obtain 


= — /no / 1 ^ a ' 4 n(n + 1 ) (4n _ 1 ) r 2 G (°>‘ 

32 

< 2 ) = — ^a (1 )2 a -4 n(n + 1 ) ( 2 n + 1 ) r 2 G<°> 
a 32 

< 


( 6 . 8 ) 


X< 2 > = ^/u.a< 1 ) a -4 r 2 G^ 0 ) 1 [n(n + l)(2n + l) , (l, 1,1, ••*,!)- 3G(° >•((), 1 , 4 , • • • , n 2 )] 
8 

X^ 2 ) = ^ v r Jl aX 1 ^ a -5/2 (n + 1) tG^ 0 ) jl + iyua -3 n(2n + 1) r 2 J • (0, 1, 2, • • * , n) . 


This formula for A.< x > is for the case k = n , i.e., for the maximization of the 
error occurring at the end of the run. Formulae (6.7) and (6.8) show that the 
coincidence of the least squares computed positions with the observed positions 
(for maximum error) that held in the first order theory (cf. (3.27) ) does not 
carry over, in general, to the second order. And, finally, the maximum squared 
error at the end of the run (k = n) for a circular orbit is 


W 2 =G(°> 


(n + 1 ) 


(> 


+ 



n 2 r 2 


s 2 


- ) 3 a" 4 n(n 4 1) r 2 ( 16 n - 1) - ^/xa -3 n 2 (n - 1) r 2 8 3 + . . . 


t 3 2 n + 1 f2 ]-i J 

(1 +|fj) 

8 n" "J j 

\ 4 -j 


(n+1)_1/2 f n + |~~ f n) 7 ( 16 n-l) -|(n-l) f 2 C 3 +...|. 
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The ratio of the second order term to the first order term for a sweep of 
one radian (f^ =1) and large n is, approximately, 


0.287 — - — 

avn + 1 

in abso lute v alue. Since in practice c can be expected to be very much smaller 
than a/n + 1 , perhaps around l/500th its size, the indication is that the first 
order theory will be adequate in the case of a single pass of moderate length. 
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